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' We set up a general formalism for models of spontaneous wave function collapse with 

["T 1 , dynamics represented by a stochastic differential equation driven by general Gaussian noises, 

not necessarily white in time. In particular, we show that the non-Schrodinger terms of the 
equation induce the collapse of the wave function to one of the common eigenstates of the 
Oh' collapsing operators, and that the collapse occurs with the correct quantum probabilities. 

We also develop a perturbation expansion of the solution of the equation with respect to the 
parameter which sets the strength of the collapse process; such an approximation allows one 



to compute the leading order terms for the deviations of the predictions of collapse models 
with respect to those of standard quantum mechanics. This analysis shows that to leading 



order, the "imaginary noise" trick can be used for non-white Gaussian noise. 
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I. INTRODUCTION 



Models of spontaneous wave function collapse [l| provide a simple consistent resolution to the 
measurement problem of quantum mechanics [2j, and at the same time provide precise indications 
for experiments which are more likely to detect possible violations of quantum linearity [3j]: these 
include e.g. fullerene diffraction experiments, decay of supercurrents, excitation of bound atomic 
and nuclear systems and several cosmological observations. The dynamics is represented by a 
stochastic Schrodinger equation of the form: 



N N 

i=i i=i 



IVft), (i) 



where H is the standard quantum Hamiltonian of the system, Ai are a set of commuting self-adjoint 
operators to whose eigenstates the wave function is driven during the collapse process, Wit are 
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N independent standard Wiener processes denned on a probability space (CI, T, P), the average 
(A{)t = (V'tl^il^) is the standard quantum expectation of A{ and 7 is a positive constant which 
sets the strength of the collapse process. 

The several collapse models which have been so far proposed differ from each other basically 
only by the choice of the localizing operators: in GRW-type models the set {A{\ corresponds 
to the set of position operators of the constituents of the given physical system, or some function 
of them; dissipative effects can be included by taking A; to be a function of both the position and 
the momentum operator of aparticle [5(, the resulting operator being non-Hermitian; in the CSL 
model for identical particles a], the index i is replaced by the space coordinate x and A(x) bacomes 
a function of the density number operator a^(x)a(x); in energy driven models Q| there appears 
only one operator A, which is identified with the Hamiltonian H. Finally, reduction models related 
to gravitational effects 8] can also be cast in the form (pQ), as shown in [9]. 

The reduction properties of Eq. ([1]) can be easily verified by computing the time evolution of 
the variance Va(£) = {A 2 )t — {A)f, of an operator A which commutes with all the operators Ai\ 
as shown e.g. in [7j, by using standard Ito calculus rules, and by setting H = 0, one gets for the 
average value Ep[V^(t)] the following equation: 

N 



E P [Vk(t)] = Va(0) - 4 7 V / dsM ¥ [Cl Ai ( s )}, 

i=i J ° 



(2) 



with CA,Ai(t) = {(A — {A)t)(Ai — (Ai) t ))t- Since the integrand on the right hand side is a non- 
negative quantity, the above relation, when applied to any operator Aj, implies that, for large times, 
the variance (t) converges to for any realization of the noise, with the possible exception only 
of a subset of Cl of measure 0; this means that any initial state \ipo) converges asymptotically, with 
probability 1, to one of the common eigenstates of the operators A{. When H ^ and moreover 
it does not commute with the other operators A4, Eq. ([1]) induces only an approximate collapse, 
the degree of approximation depending on the relative strength of the Schrodinger term and of the 
collapse terms which define the equation. 

A very useful mathematical property of Eq. ([1]) is that its physical predictions concerning the 
outcomes of measurements are, in terms of statistical expectations, invariant under a phase change 
in the noise. As a matter of fact, let us consider the following class of stochastic Schrodinger 
equations: 
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l -Hdt + V7"5>^ " ^( A i)t)dWi,t ~ \ ^(|£| 2 4 2 " 2^MAi) t + en(A t ) 2 t )dt 

i=l i=l 



(3) 
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where £ = £ R + «£i is a constant complex factor; of course, when £ = 1 we recover our original 
collapse equation. An easy application of Ito calculus leads to the following equation for the density 
matrix p(t) = E P [\ip t ) {ip t \]: 

, ■ N 

f t p(t) = -~[H,p(t)] + l\tfY,(2A P (t)A-{Alp(t)}), (4) 

i=l 

which is of the Lindblad type and has the remarkable property that it depends only on th e sq uare 
modulus of £. Since, within collapse models, the statistics of the outcome of experiments [10] can 
be expressed by the averages Ep[(^t|O|-0t)] = Ti[p(t)0], where O is a self-adjoint operator, we see 
that in order to compute experimental predictions, one can use in place of Eq. ([T]) any stochastic 
equation of the type ([3]) which satisfies the constraint |£| = 1; in most cases, it is convenient to 
choose the equation corresponding to £ = i, since it is linear, thus much easier to solve [llj. Of 
course this does not mean that all equations of the form ([3]), having the same value of |£|, are 
equivalent; on the contrary, in general they generate completely different evolutions for the wave 
function. For example, when £ = 1, as we have seen, the corresponding equation induces collapse 
of the wave function, since the property (|2|) for the variance holds true; while on the other hand 
with £ = i, the corresponding equation is linear, thus the wave function does not collapse in this 
case. This notwithstanding, all quantities of the form Ep[(V^|0|V^)] turn out to be the same for 
the two equations, and for similar ones with complex £ of modulus unity. 

The aim of this paper is to generalize Eq. (pQ) in order to include also non-white, Gaussian, 



13, Q, fl, but a 



stochastic processes. Some results have already appeared in the literature 
general analysis is still lacking, mainly because the new dynamics is not Markovian and thus is 
more difficult to describe mathematically. There are two main reasons why one should consider 
collapse models driven by non-white noises. First of all, it is important to understand how the 
collapse mechanism and other physical properties, such as the time evolution of the mean energy, 
depend on the type of noise driving the collapse of the wave function. Since these properties 
are directly connected to physical predictions which differ from those given by standard quantum 
mechanics, such differences, and thus the possibility of experimentally testing collapse models, 
could significantly change depending on the type of noise entering the collapse equation. The 
second reason for such an analysis is that a non-white noise, unlike a Wiener process, can be 
identified with a physical field; accordingly, one can try to connect the collapse mechanism to some 
other physical process occurring in Nature, possibly having a cosmological origin; we will come 
back to this point in the final section. 

This paper contains two main sections, which set up the general formalism for non-Markovian 



collapse models by following two different paths: in Sec. [II] we follow the same argument used in [f| 
to derive, from a generic diffusion process in Hilbert space, Eq. ([I]) as the correct collapse equation; 
in Sec. IIIII instead we follow the strategy used in 16| to obtain the same equation from general 



requirements on the dynamics of the density matrix. We will prove the two following main results: 

1. We will show that, if one neglects the quantum Hamiltonian H, the dynamics leads to the 
collapse of the wave function to one of the common eigenstates of the localizing operators, 
with the correct quantum probabilities. 

2. We will develop a perturbation expansion of the solution of the equation with respect to the 
coupling constant ^7; when applied to experimental predictions on microscopic systems, it 
provides the leading order term for the deviation of such predictions from those given by 
standard quantum mechanics. 

Concerning this second result, we will also show that, at least to order 7, the equation for the 
statistical operator depends only on the absolute value of £, precisely as discussed before for the 
white-noise case; this means that, to this order, one can employ the useful trick of replacing the real- 
noise (£ = 1), non-linear collapsing equation with an imaginary-noise (£ = i), linear non-collapsing 
equation, thus considerably simplifying calculations. 

Throughout the calculations, one has to make sure that, at any stage, one recovers the correct 
white- noise limit; however, one has to keep in mind that the white-noise limit of the non- white 
collapsing equation will not be given by (pQ) or ([3]), since these are ltd equations; instead, as 
explained e.g. in [13], an equation containing non-white noises reduces, in the white- noise limit, 
to a Stratonovich equation. Accordingly, we write the Stratonovich equation corresponding to the 
ltd equation ([3]), which is Q]: 

N 



d\il> t 
dt 



i=l 
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(5) 



In the second expression, we have written the third term in square brackets in a form that corre- 
sponds to the white noise limit of Eq. (fT8|) below, with the remainder of the order 7 part playing 
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the role of a shift in the mean value of the noise, corresponding to a change of measure, as dis- 
cussed in Sec. IV below. Of course, whether one uses Eq. ((3]) or ([5]) , the corresponding equation 
for the statistical operator p(t), which determines the evolution of statistical ensembles of states, 
will always be given by (JH). 



II. LINEAR AND NONLINEAR COLLAPSE EQUATION 

Following the path outlined in [6] to construct the continuous generalization of the GRW model 
in terms of an Ito stochastic differential equation of the type (pQ), let us consider a diffusion process 
for the wave function in Hilbert space having the form: 

dm)) 



N 

~h + V7C^A^(t) + o 

" i=l 



\m, (6) 



dt 

where, as before, H is the standard quantum Hamiltonian of the system, Ai are commuting self- 
adjoint operators, 7 is a positive coupling constant, £ = £ R + i6 is a constant complex factor, while 
O is a linear operator yet to be defined 1 . The noises Wi(t) are real Gaussian random processes 
defined on a probability space (O, J 7 , Q) whose mean and correlation functions are, respectively: 

E Q [ Wi (t)} = 0, Eqlwiit^Wjih)) = D ij (t 1 ,t 2 ). (7) 

When £ R 7^ 0, which is the case for collapse models, Eq. © does not preserve the norm of the wave 
function; therefore we introduce the normalized vector: 



(assuming of course that the norm of \4>(t)) does not vanish) which we take as the physical vector 
describing the random state of the system at time t. 

The measure Q previously introduced is not the correct physical probability since it does not 
lead to a collapse that respects the Born probability rule; the right physical probability, which we 
shall call P, is defined as follows: 

P[F] = E Q [l F (0(tMt)>] VFef, (9) 

where lp is the indicator function associated to the measurable subset F of fi. This definition 
corresponds to the assumption of the GRW model, according to which a collapse (in space) is more 
likely to occur where the wave function is larger, as postulated by the Born probability rule. 



1 As we shall see, the operator O will not be a standard linear operator, since it will act on a vector also through 
its dependence on the noises Wi(t), by means of functional derivatives. In this respect, our approach differs from 
the standard one based on Ito's diffusion equations in Hilbert spaces. 
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To summarize, an initial state |^>(0)) is driven by the stochastic dynamics into an ensemble 
of states \ip(t)) of the form ([8]), where \4>(t)) solves Eq. ([6]); the distribution of states within the 
ensemble is given by the probability P defined in ([9]). 

Note that the definition ([9]) of P not only matches with the Born probability rule, but is also 
necessary in order to prevent the possibility of usin g th e collapse mechanism to send information at 
a speed faster than the speed of light. As shown in [Ha], when the dynamics of a statistical operator 
p(t) is nonlinear, it is in general possible to send faster-than-light signals which can be used by two 
spacedike separated observers to communicate with each other; this possibility is instead forbidden 
when the evolution is linear. In our case, according to the previous assumptions, p(t) is defined as 
follows: 

p(t) = E P [|^(t))^(t)|]; (10) 

but according to Eqs. ([8]) and ((9j) we have the mathematical equality: 

E P [|^(t))(^(t)|] = EQ[|0(t))<0(t)|], (11) 

and since \(j>(t)) solves the linear equation ©, it follows that the operator which maps p(t) to 
p(t + dt) is linear, thus not allowing faster-thandight signalling of the type just discussed. 

We still need to prove that Eq. ([9]) correctly defines a probability measure; it is easy to check 
that all properties are satisfied expect for the normalization condition P[f2] = Eq[(<^|<^)] = 1, 
which in general is not fulfilled unless the operator O takes a particular form. In fact, by using the 
Furutsu-Novikov formula |ig| ] : 

V r SF[{w(t)}Y 



EQ[F[{w(t)}]wi(t)] 



V / dsD ij (t,s)E i 
,=i Jo 



5vjj(s) 



(12) 



,r- 

which holds for a generic functional F[{u?(i)}] of the Gaussian noises Wi(t), i = 1, . . . N satisfy- 
ing ([7]), and computed from initial time to time t, one can immediately prove that: 

, N ft jC 

-EqK&I^)] =0 if: O = -2V7£rE^ / dsD l3 (t,s)—— (13) 
dt ^ J 5 Wj {s) 

Accordingly, the linear equation which, together with Eq. @, induces collapse of the wave function 
with the correct quantum probabilities and, at the same time, does not allow one to use the collapse 
process to send signals at faster-thandight speed, is [l^ : 

N N 



dim 

dt 



■~H + V7£ Y^Awiit) - 2^7 £r Y. A i[ dsDij&s) 
i=i i,j=l Jo 



dlVj(s) 



W))- (14) 
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As foreseen, this equation is non-Markovian and for this reason is highly non-trivial, since the 
future evolution, which involves the whole past, depends on the combined effect of the standard 
Hamiltonian H and the collapsing operators Ai. This dynamics is not easy to unfold if these 
operators do not commute among themselves, as is usually the case. 

The equation for the normalized vector \4>(t)) does not have a closed form, unless the functional 
derivative of \4>(t)} can be explicitly computed; as we shall se in the next sections, this happens 
when the Hamiltonian H is neglected (or when it commutes with the operators Ai), and when one 
writes the evolution as a perturbation expansion with respect to the relevant parameters. 

Before moving on, we make a few comments about the change of measure as defined in ([9]). 
In the white noise case one can prove, under suitable hypotheses on the operators H and Ai 
(see e.g. that (^>(t)\<p(t)) is martingale with Q-mean equal to 1; this ensures that one can 

consistently use (<f)(t)\<f)(t)) as a Radon-Nikodym derivative of a new probability measure P with 
respect to Q, as we have assumed more heuristically in the previous paragraphs. Here we will not 
attempt to prove that {(fi(t)\(f)(t)) satisfies the required properties also in the more general case 
of non- white Gaussian processes, leaving the analysis of the conditions under which this is true 
to future research. Secondly, Girsanov's theorem provides, in the white noise case, a connection 
between Wiener processes with respect to the measure Q and Wiener processes with respect to the 
transformed measure P. It would be interesting to see whether a similar theorem can be proved also 
in the non-white noise case, and whether Q-Gaussian processes can be connected to P-Gaussian 
processes; we will come back on this point in Sec. IV. 



A. Collapse of the state vector 



We now show that, when the standard quantum Hamiltonian H is set to 0, the dynamics 
induces the collapse of the state vector \ift(t)) to one of the common eigenstates of the operators 
Ai. As shown in [14j |. if H is neglected so that all operators entering the equation commute among 
themselves, then the functional derivative can be explicitly computed and Eq. (I14B reduces to: 



dim 

dt 



N N 

^ii^AiW^t) - 2 7 ££r AiAjFi.it) 

i=l i,j = l 



(15) 



where we have defined: 



Fij{t) = [ dsD i3 \t,s). 
Jo 



(16) 
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That Eq. ()15[) is equivalent to Eq. (|14p in the limit H = can easily be seen by integration of 
Eq. (I15p . which is trivial since all operators commute, from which one obtains the relation: 

ft^) W)> = Vt^M*)), s<«, (17) 

a relation which will often be used in the following calculations. 

The equation for the normalized vector \ip{t)) can now be directly computed from the defini- 
tion (ED: 



d\m) 

dt 



N N 



W)h (18) 



with the expectations (...)t computed in the state \ip(t)), that is (O)t = (V ; (i)|0|V'(^))- 

The equation for the statistical operator p{t) can now be computed either from Eq. (|18|) through 
the definition of Eq. (jlpp or from Eq. (|15p through the equivalence of Eq. (jlip : in both cases one 
gets: 

^ p(t) = 7 l^ 2 J] (^ A i + A i A i ~ A ^ A iP^) ~ p(t) A i A i) Fij{t), (19) 

which correctly reduces to ([4]) in the white-noise limit. 2 

We have now all the necessary formulas to compute the time evolution of quantities such as 
Ep[(^4 n )j] and Ep[(A)™], where A is a self-adjoint operator commuting with all the operators Ai. 
Because of the relation Ep[(-A n )t] = Tr[A n p(t)], which is a consequence of Eq. (flOl) . and because of 
the trace-preserving structure of Eq. (|19p . one immediately has: 

j t m( An )t] = o. (20) 

By using the change of measure and through a direct calculation, one finds for Ep[(A)™] instead 
that: 

j t Mi A )t\ = 2^Um( A )r 1 (n(AX) t -(n-l)(A) t (X) t )], 

N N 

x = y, A i w i(t) - 2 Vi ?r A i A 3 F ij(t). (21) 

t=l i,j=l 



2 Note that in the white noise limit, through Eq. (|16[) one encounters the integral of a delta function at the endpoint 
of an interval, which is L dsS(t — s) = 1/2, since 5(t) = S(—t) and dtS(t) = 1. This enters both in comparing 
the white noise limit of Eq. (Tl9|) to Eq. (j4|, and the white noise limit of Eq. (fl8|l to Eq. ([5]). 
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The terms proportional to Wi(t) can be rewritten by using the Furutsu-Novikov formula, together 
with the equality: 



y/l[e{AjO) t + Z{OAj) t - 2UO)t(A J } t ], (22) 



Swj(s) 

which is valid for any operator O, and which can be directly proved from the definition (O)t = 
(4>(t)\0\(t)(t)) / {(j)(t)\(f)(t)} , together with Eq. (fTTl) . After a rather lengthy calculation, one can prove 
that Eq. (|2ip simplifies to: 



i-3 

We now apply Eqs. (f20j) and ([23]) to compute the time evolution of the variance Vk(£) = 
{A 2 ) t — (A) 2 of the operator A; we obtain: 

Ep[Vk(*)] = V A (0) - 8^7^ / dsE F [((A i -(A i ) s )A) s ((A j -(A j ) s )A) s ]F ij (s). (24) 

i.j 

Now, the same argument used in Eq. to prove the reduction of Eq. ((TJ) holds true: when the 
matrix Fij(s) is positive definite in the limit as t —* oo, Eq. (|24|) is consistent if only if, for large 
times, ((A; — (Ai)t)A) t goes to zero for any i almost surely (a.s.) (i.e., except on a subset of 
of realizations of the noise of P- measure 0); in particular, if we take A equal to any one of the 
operators Ai we have: 

hm [<^) 4 _ ( At >2] = hm V Ai (t) = a.s. V i, (25) 

t— >oo t— >oo 

which is the desired result. Moreover, due to Eq. (|20p . the average value of {P ai )t remains constant 
in time, with P ai the projector on any eigenspace of Ai with eigenvalue dj, which means that the 
collapse occurs with the correct quantum probabilities. 



B. Perturbation expansion to order 7 

The approximation used in the previous subsection, which consisted in neglecting the quantum 
Hamiltonian H, is useful when the system under study is macroscopic, since in this case the effect 
of the collapsing terms is typically much stronger than that of H: this is precisely the reason 
why collapse models ensure the localization of the wave function at the macroscopic level. For 
microscopic systems, on the other hand, such an approximation is no longer valid; just the reverse 
from the macroscopic case, at the microscopic level the effect of the collapsing terms represents 
typically only a small perturbation on the standard quantum evolution: for this reason collapse 
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models agree very well with standard quantum mechanical predictions. It then becomes meaningful, 
for micro-systems, to perform a perturbation expansion of the evolution of the state vector with 
respect to the parameter yfy, in order to compute the leading terms representing the deviations of 
the predictions of collapse models from those given by standard quantum mechanics. To this end, 
let us introduce the interaction picture operators and states: 



Mt) = U\t)AiU{t) 
Eq. (|14p then becomes: 

d\m) 



9(t)) = U\t)\m, U(t) = exp(~Ht); (26) 



dt 



N N t 

^ Y^Ai{t) Wi (t) - 2^7 & Y\ Mt) / dsD^s 

i=l Li=l J ° 



Swj (s) 



m- (27) 



The perturbation expansion of |<^ ! (i)) with respect to the parameter ypy reads: 

w{t)) = + v^WM) + 7i0 i 2 (t)) + ..., 

while the functional derivative acts on \4> l {t)) as follows: 

d :l^(*)> = >/7t4tI^(*)) + iT^m)) + ---; 



(28) 



(29) 



5vjj(s)' T w/ v ' Swj(s) 8wj(s) 

the second term must be of order 7, since the functional derivative brings down a term proportional 
to yp); for the same reason, the third term is of order 7 3 / 2 . This means that the perturbation 
expansion can be explicitly carried out, despite the functional derivative appearing in Eq. (f27|) . 
and we get the following results to order 7: 



order 0: ^|#>(*)> = 0, 

d N 
order ^7: — = ^(*K(*)l0o), 

i=l 

AT 



(30) 
(31) 



d N N rt 

order 7: jM(t)} = ^ M^M^M®) ~ 2 ^ Yl / dsj4i ( s ) D ij (t,s)\ . (32) 

i=l i,7=l 



Going back to the Schrodinger picture, Eq. (fTi|) . to order 7, reduces to: 

N 



d_ 
~dt 



mm + 



i=l 



+ 7^5" dsA i A j (s-t)(Zw i (t)w j (s)-2£ R D ij (t,s)) 
rr* 1 ^ Jo 



(33) 
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with \4>(t)} the total wave function, and \<f>o(t)) its zeroth order part. By use of Eqs. (30)-(32), this 
can also be written entirely in terms of \4>(t)) as 

'o 



j:H + ^7 AiWi ^ ~ 2 7^hX /„ dsAiAjis- t)D i3 {t, s) 



i=l 



(34) 



We can now see more clearly the non-Mar kovian nature of the evolution, because of the presence 
of the term Aj{s — t), which depends on the past effect of H on Aj. When H = 0, one immediately 
sees that by using Eq. (fTBj) . the \4>{t)) evolution equation of Eq. (f34l) reduces to the evolution 
equation given in Eq. (|15p . 

The corresponding equation for the normalized vector \ip(t)) defined by Eq. ([5]) can now be 
obtained by a straightforward application of the chain rule for differentiation, with the result 



N 



-(H + ih 7 AS A) +y/l^2^Ai - UAi)t)Wi(t) + j(0 S A ~ (0 SA ) 



i=l 



\m), (35) 



again with the expectations (...)t computed in the state \i/j(t)}. Here Oasa and Osa are respectively 
the anti-self-adjoint and self-adjoint parts of the operator O defined by the final term inside the 
square brackets of Eq. (|34p . and are given explicitly by 

N 



rt 

Oasa = -Y] / - t)j + ^? a {ii, - i)}) s) (36) 

i,i=i Jo 

N r t 

S A = -Y] / ds(e R {A l ,A j (s - t)} + iCrUAi,Ms - t)])Z>y(t, S ). (37) 
id=i Jo 



The equation for the statistical operator can now be computed by resorting to relations (|10p 
and (jlip . The calculation proceeds most directly from the \4>(t)) evolution equation given in 
Eq. (|33|) . since in this equation all dependence on the noise is explicit. To the order to which we 
are working, we can replace the 7 term in this equation by its expectation Eq, giving the simplified 
evolution equation 

v N rt 

'0 



± t \ m = - i - Hm) + 



Vl^AiWiit)-^] 2 ^ IJsAiAjis-^Dijfas) 

M=l 



i=l 



*o(t)>. (38) 



It is then straightforward to compute EQ[|^>(t))(0(i)|], with the result 3 



3 P. Pearle has pointed out that, in the real noise (£ = 1) case, this equation follows from differentiation of the 
integrated expression given in Eq. (4.11) of his Physical Review article cited in [l^], which he suggests is exact 
when time-ordering is included. 
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-AiAj {s - t)p(t) - p(t)Aj (s - t)Ai] . (39) 

As we can see, the above equation, which is correct to order 7, depends only on the absolute value 
of £; this means that all physical predictions are, at least to order 7, independent of the phase of 
£, precisely as in the white-noise case. As a consequence, in order to compute physical predictions, 
one can resort to Eq. (|14p with £ = i, which is much simpler since the integro-differential term 
vanishes, and one is left with a standard Schrodinger equation with a random hermitian potential. 



III. AN ALTERNATIVE CONSTRUCTION OF THE NONLINEAR COLLAPSE 

EQUATION 

We give in this section an alternative construction of the nonlinear collapse equation in the 
case of non-white Gaussian noise. Instead of starting from a linear equation and using a change 
of measure, we work throughout with a norm-preserving nonlinear equation, and use perturbation 
theory to determine the structure of the order 7 term which guarantees that the corresponding 
density matrix evolution will have the Lindblad form [21] in the Markovian limit. 
Thus, returning to Eq. ([6]), we now start from 

N 

h 



dim) 



-~H + ^^ A i - U(Ai)t)Mt) + l{B SA - (B SA ) t ) + jB ASA 
i=l 



\m), (40) 



dt 

with B$a an d Basa respectively a self-adjoint operator and an anti-self-adjoint operator to be 
determined, and with Wi{t) now a non- white Gaussian noise obeying 

%K(i)] = 0, E M MtiH-(i 2 )] = DyQuh) (41) 

with respect to the measure R. We shall determine B$a and Basa to simplify the evolution 
equation for p(t) = Er [/}(£)] = Ejj[|'(/'(t))(^(t)|], in such a way that in the Markovian limit it 
reduces to a Lindblad evolution in the first standard form, 

<W) r /.x 

= - ~[H,p{t)\ +Y,ai3 - \ {F]Fi,P(t)}) > (42) 

with Fi suitable functions of {Aj}, and with the coefficients oy determined by the noise expectation 
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By construction, Eq. (|40p preserves the normalization of the state vector \ip(t)) under time 
evolution. Prom this equation, and its adjoint, one easily finds that the pure state density matrix 
p(t) obeys the evolution equation 

d ^ L =- kH,p(t)] + [^i^^Mt)+lB ASA ,p(t)] 

i=l 

N 

+{V7£r E^ " {Ai)t)wi(t) + 1 {BsA ~ {B SA ) t ),p{t)}. (43) 
i=i 

Taking the expectation of this, and retaining terms through order 7 but dropping terms of order 
7 3 / 2 , we get 

N N 
<u UH,p(t)] + V7i&£[A,%[£(tH(*)]]+7^ 

i=l i=l 



dp{t) _ i 

2 V7£r Yl Mp(t)(Ai) t Wi(t)] + 7 {B SA ~ (B SA )t, p(t)} ■ 



N 



i=l 



We now use the Furutsu-Novikov formula 



E R [F[{w(t)}} Wl (t)} = V f dsD lJ (t, S )E R \ 6 -^j^ ] 

~[ Jo L dw j { s > 



(44) 



(45) 



first with F = p(t), and then with F = (Ai)tp(t) = p(t)Trp(t)Ai. By the chain rule, 

dwj{s) owj{s) 0Wj{s) 

so for both choices of F in Eq. ([15]) what we need is 5p(t)/Swj(s), calculated through terms of 
order ypy. This can be calculated directly by integrating the differential equation of Eq. (|40p . In 
terms of the interaction picture operators 

Aj(s -t) = ei^-^Aje-i 11 ^, (47) 



a simple calculation gives 
Sffl 



VlHAMs - t),p(t)] + UiMs - t) - {Aj)„ p(t)}} . (48) 



Swj(s) 

Substituting this into Eq. (|44p gives a lengthy expression, which on algebraic simplification, and 
noting that (Aj(s—t))t = (Aj) s , gives a result that may be summarized as follows. Let us abbreviate 

N rt 

(49) 



S ij = ^ / dsDij(t, s), 
i,i=l Jo 
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so that Sij acting on a function of t, s gives a function only of t. Then we find 
dp(t) i 



dt h 

-2 , t 2 



[H,p(t)] 

7(£r + &)Sij[Aip{t)Aj{s -t) + Aj(s - t)p(t)A l - A l A j {s - t)p(t) - p{t)Aj{s - t)A l ] 
+l[Si j & lj + ittUCij) + Basa, p(t)] 

+l{SiM(. C H ~ ( C v)t) + iUn(Fij ~ (Fij)t)} + B S A ~ (B S a)u p(f)}, (50) 
where we have introduced the condensed notations 
Fij =[Ai,Aj(s-t)], 

Cij ={^,^-(s - t)} - 2A i (A j (s - t)) t - 2(A i ) t A j (s - t). (51) 
Let us now make the following choice of the previously undetermined operators B$a and Basa, 

BsA = - Sij{£,l.Cij + itii^nFij), 

Basa = - S lJ (ClF, + i^dj). (52) 

Then the final two lines of Eq. (|50p cancel to zero, and we are left with 
dp(t) i N rt 



~[H,p(t)]+7(Zl + g) E / dsD tJ (t,s) 

i,j=l Jo 

x [Aip(t)Aj (s-t)+ Aj (s - t)p{t)Ai - AiAj (s - t)p{t) - p(t)Aj (s - t)Ai] . (53) 



Thus, we recover from this approach the evolution of Eq. (|39p of Sec. II. In the Markovian limit in 
which D{j(s, t) decays rapidly when s is not close to t, we can approximate s — t ~ in Eq. f)53[) . 
We then have Aj(s-t) ~ Aj, which by assumption commutes with Ai, and the evolution of Eq. (|53p 
takes the first standard form of a Lindblad evolution, 



dp(t) i 



N ft 

^[#,p(t)]+ 7 (£ + £?)E / dsDij&s) 
xlAip^Aj + A jP {t)Ai - {AiAj,p(t)}\. (54) 



Eq. (|54p also follows, without the Markovian approximation, if one neglects the Hamiltonian H, 
since then the interaction and Schrodinger pictures coincide, and one has exactly Aj{s — t) = Aj. 
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IV. CHANGE OF MEASURE AND OTHER CONCLUDING REMARKS 

A. Change of Measure 

In Sees. II and III we have given two different derivations of the time evolution equation for 
the normalized wave function \ip(t)}, and of the corresponding evolution equation for the density 
matrix p(t). By construction, the density matrix evolutions of Eq. (|39p and Eq. (|53p are the same. 
However, a comparison of the operators Osa and Oasa of Sec. II with the operators B$a and Basa 
of Sec. Ill shows that these are not the same, and hence the corresponding evolution equations 
for \ip{t)) of Eq. (|35p and Eq. (|40p are not the same. This means that our two constructions 
correspond to inequivalent unravelings of the same density matrix evolution: in general the noises, 
expectations, and wave functions \ip(t)) of Sec. II and Sec. Ill are not the same, even though they 
lead to the same evolution equation for the noise-averaged density matrix. 

However, in certain special cases the formulations of Sees. II and III are related by a time- 
dependent shift in the mean values of the Gaussian noise variables. In order for the functions 
\ip(t)) used in the two derivations to be identical, they must obey the same time evolution equation. 
Changing notation, by denoting the noise of Sec. Ill as Wi(t), we find that the evolution equations 
for \ip(t)) of Eq. ([35]) and Eq. ([40]) become identical when the noises Wi(t) and Wi(t) are related to 
leading order in y/j by 

N N JV _t 

^Wi(t)Ai = ^2wi(t)Ai - 2^7 £r Yl / dsDijfas^AiiAjis - t)) t + {A^A^s - t)}, (55) 

1=1 i=l i,j=l ^ 

with 

®&[wi(t)M s )] = EqK(*H-(s)] = Dij(t, s) (56) 

to zeroth order in ^Py. 

Equation (|55|) is consistent only if the operators Aj(s — t) can be expanded over a basis of the 
time zero operators Aj, with c- number coefficients, that is 

N 

A J (s-t) = Y,K ji (s-t)A i . (57) 
i=i 

This is automatically true (i) when the Hamiltonian H vanishes, since then Aj(s — t) is time 
independent, and (ii) in the white noise limit, since then Dij(t, s) oc 6(t — s), and so Eq. (|55p only 
involves A/(0) = Aj. It is approximately true (iii) whenever a Markovian approximation to the 
time evolution is valid. In general, however, Eq. (|55l) involves on the right hand side operators 
Aj(s — t) that are linearly independent of the operators Aj, and so cannot be satisfied. 
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When Eq. (|57p holds, then Eq. (|55p simplifies to take the form 

N N N 



J2wi(t)Ai = ^Wi{t)Ai -YjKi(t)A h (58) 

i=l i=l i=l 

with -ftTj(t) given by 

N t N 
Ki(t) =2^7 & / ds tAi(*, a)<^(a - t)) t + D im (s - t^A^K^s - t)\. (59) 
j=l ^° m=l 

Eq. (|58p can clearly be satisfied by making a c- number shift in the noise variable for each i, 

Wi{t) = Wi (t) - Ki(t) (60) 

The relation at time t between the measure R of Sec. Ill, in which Wi(t) has zero mean, and the 
measure Q of Sec. II, in which Wi(t) has zero mean, is then to first order in y^y, for any argument 
O, 

E M [0]=M Q [W(t)0]. (61) 
Here the weighting factor W(t) is given by 

(62) 



W(t) = 1 + V / daCi(t,8)wi(8) 
i=l J o 



where Ci(t,s) obeys the integral equation 

N rl 



K i( t )=Y J l dsD i:j (t, 8)Cj{t,8) 

i=l J ° 



(63) 



This choice of C{(t, s) guarantees that when O is taken as Wi in Eq. (|6ip . one finds that Ejgfwj] = 0, 
as needed. An analogous c-number shift of the noise variable enters in comparing Eq. ([5]) of Sec. 
I with the white noise limit of Eq. (|18p in Sec. IIA. 



B. Final Remarks 



In the preceding sections we have shown that the white noise formalism is robust under a 
generalization to the physically more realistic assumption of non-white noise. Both the proof of 
state vector reduction, and the "imaginary time" trick for calculating physical effects of the noise, 
carry over to the non-white noise case, to leading quadratic order in the noise strength. 
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We wish here to elaborate on some implications of our calculations for models of the non- 
white noise. We recall that the noise autocorrelation Dij(t\,t2) is defined as the expectation 
^Q[wi(ti)wj(t2)], and that the condition for state vector reduction is that the time integral 



F i:j (t) = [ dsDijfas) 
Jo 



(64) 



should be a positive definite matrix in the limit as t — > oo. 4 Let us now assume time translation 
invariance, which implies that Dij(t, s) = Dij(t — s), and investigate what the requirement for state 
vector reduction means in terms of the spectral decomposition of D^j. Writing 



Dij(t 



duj^ij{oj) cos u(t 



we have 



sin tot 



00 du 

— sin wjij(u/t). 



10 w jo u 

Thus, assuming that j(uj) is smooth in the neighborhood of u> = 0, we find 



f°° du 

lim Fij(t) = 7i,(0) / — sinu = 7i,-(0)7r. 
t->°° Jo u 



(65) 



(66) 



(67) 



Hence, the reduction requirement is satisfied when 7^(0) is a positive definite matrix in In 
particular, the spectral weight 7ij(w) can have a cutoff at a finite upper limit lu = CL> max , without 
in any way affecting the reduction argument. The possibility of such an upper cutoff has been 
discussed in the review of Bassi and Ghirardi [l| , and as noted by Adler and Ramazanoglu 22|] , 



is suggested on physical grounds by existing upper limits on noise-induced gamma ray emission. 

Let us next consider the rate of secular energy increase induced by the noise, taking advantage of 
the "imaginary noise" trick to write the noise term as an addition to the Hamiltonian. We consider 
the simple model describing a particle of mass m moving in one dimension, with a non-white noise 
coupling to its coordinate x, with Hamiltonian 

„2 



H 



EL 

2m 



Cwtx. 



(68) 



Because the noise term does not commute with p, the kinetic energy p 2 /(2m) increases over time. 
We have for the expected rate of energy gain, 



dt 



- p2 - 


= m _1 E 


dp 






2m 





(69) 



4 More generally, as can be seen from Eq. (|24p . the reduction requirement is satisfied when J Q dsFij(s) — » 00 
t — ► 00. 
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From the Heisenberg equations of motion implied by Eq. (|68p we find dx/dt = p, dp/dt = Cwt, 
and so 



p(t) = p(0) + C I duw. 



Hence the rate of energy gain is 
d 



dt 



-E 



\ p2 ~ 


_ c 2 r 


2m 


™ Jo 



t q2 

duK [w t w u ] 



m 



t (ji 

duD(t,u) = — F(t), 
m 



(70) 



(71) 



and therefore is governed by the same integral over the autocorrelation function as the reduction 
rate. Thus, when the reduction condition lim^oo F(t) > is obeyed, the rate of noise-induced 
energy production is nonzero at large times, a result which readily generalizes to more realistic 
reduction models. As reviewed in 3|, this leads to various upper bounds on the noise strength. 
This conclusion can be evaded in the generic case of multiple operators Ai, if the only nonvanishing 
eigenvalue of Fij(t) as t — > oo is the one associated with the total energy. An interesting model 
where this case is realized, but in which localizing reduction still occurs in an approximate sense, 
is given by taking Dij{t, s) to be a correlation function associated with thermal noise, as might be 
expected if state vector reduction is induced by some type of cosmological relic field. A detailed 
examination of the thermal noise model will be given elsewhere 23]. 
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